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GAUSSIAN PROCESSES FROM A UNIQUE POISSON PROCESS 
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Abstract. We consider two independent Gaussian processes that admit a 
representation in terms of a stochastic integral of a deterministic kernel with 
respect to a standard Wiener process. In this paper we construct two families 
of processes, from a unique Poisson process, the finite dimensional distributions 
of which converge in law towards the finite dimensional distributions of the two 
independent Gaussian processes. 

As an application of this result we obtain families of processes that con- 
verge in law towards fractional Brownian motion and sub-fractional Brownian 
motion. 



1. Introduction and preliminaries 

Let f{t, ■) and g{t, •) be functions of L'^(R+) for all te [0, T], T > and consider 
the processes given by 

(1.1) Yf = y^°" f{t,s)dWs,te[0,T]^ 
and 

(1.2) ^'"{i 9{t,s)dWs,te[0,T] 

where W — {Ws, s > 0} and W — {Ws, s > 0} are independent standard Brownian 
motions. 

The aim of this paper is to construct two families of processes, from a unique 
Poisson process, that converge, in the sense of the finite dimensional distributions, 
to the processes and Y^. We will use this result later in order to prove weak 
convergence results towards different kinds of processes such as fractional Brownian 
motion and sub-fractional Brownian motion. 

It is well known the result by Stroock (see |Str82j ) where it is shown that the 
family of processes 



N s 

(-1) 7^ ds, t e [o,r] 



1 N s 

defined from the kernels 0^ = -(—1) ^ , converges in law in C([0, T]) to a standard 
Brownian motion, where N = {Ns, s > 0} is a standard Poisson process. This kind 
of processes were introduced by Kac in [Kac74] in order to write the solution of 
telegrapher's equation in terms of Poisson process. 

On the other hand, Delgado and Jolis (see |DJOO| ) extend this result to processes 
represented by a stochastic integral, with respect to a standard Wiener process, of 
a deterministic kernel that satisfies some regularity conditions. 
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A generalization of Stroock's result can be found in j Bar01| . where it is proved 
that the family 

(1.3) I^xt{t) = -J\''''iids, ie[0,T]| 

converges in law in C([0, T]) to a complex Brownian motion, for 9 G (0, tt) U (tt, 27r). 
Particularly, the real part and the imaginary part of (|1.3p tend to independent 
standard Brownian motions. 

In this paper, given {Ns,s > 0} a standard Poisson process and 9 G (0, tt) U 
(tt, 2tt), we consider the following families of approximating processes 

(1.4) Y/^|J^ /(i,s)cos(0iVa.)ds, tG[0,T]| 
and 

(1.5) ^"'"{ei g{t, s) sin {9Nz^)ds, tG[0,r]|. 

The main result of this paper is the proof that the finite dimensional distributions 
of the processes Y/ and V/ converge in law to the finite dimensional distributions 
of the processes and Y^ given by (jl.ip and (|1.2p . respectively. 

It is important to note that the processes Yj and Yf are both functionally 
dependent. Nevertheless, integrating and taking limits, we obtain two independent 
processes. 

As an application of this result it can be obtained approximations for different 
examples of centered Gaussian processes, among others, fractional Brownian motion 
and sub-fractional Brownian motion. 

Recall that fractional Brownian motion (fBm for short) = {B^ (t), i > 0} is 
a centered Gaussian process with covariancc function 

(1.6) Cov(i?f,i?f)^l(s«+t«-|.-r) 

where H e (0, 2). Usually fBm is defined with Hurst parameter belonging to the 
interval (0, 1) with the corresponding covariance, but in order to compare it with 
sub- fBm we use the stated representation with H G (0, 2). 

On the other hand. Sub-fractional Brownian motion (sub-fBm for brevity) = 
{S^{t), i > 0} is a centered Gaussian process with covariance function 

(1.7) Cov(5f , 5f ) ^s" + t"-^[{s + t)" + \s~ t\"] 
where H G (0,2). 

This process was introduced by Bojdecki et al. in 2004 (see jBGT04j ) as an 
intermediate process between standard Brownian motion and fractional Brownian 
motion. Note that both fBm and sub-fBm are standard Brownian motions for 
H = \. 

For H ^ I, sub-fBm preserves some of the main properties of fBm, such as long- 
range dependence, but its increments are not stationary; they are more weakly 
correlated on non-overlapping intervals than fBm ones, and their covariance decays 
polynomially at a higher rate as the distance between the intervals tends to infinity. 
For a more detailed discussion of sub-fBm and its properties we refer the reader to 
[BGT04j . Some properties of this process have also been studied in |Tud08| and 
[Tud07| . On the other hand there is an extension of sub-fBm in [BGT07) . 
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In | RdCT09] (see Theorem [33] below) the authors obtain a decomposition of the 
sub-fBm in terms of fBm and another process with absolutely continuous trajecto- 
ries, X" = {X^,t > 0}, which is defined by Lei and Nualart in |LN09| by 

(1.8) ^f=/ (l-e-'^*)r-^dW, 

where is a standard Brownian motion. Lei and Nualart introduce this process 
in order to obtain a decomposition of bifractional Brownian motion into the sum 
of a transformation of and a fBm. 

The decomposition is different for H £ (0,1) and H S (1,2). In the first case, 
sub-fBm is obtained as a sum of two independent processes, a fBm and the process 
defined by (|1.8p . while for H G (1,2) is fBm that is decomposed into the sum of 
the process (11.81) and a sub-fBm, being these independents. 

The paper is organized as follows. In Section 2 we will prove the general result 
of weak convergence, in the sense of the finite dimensional distributions, towards 
integrals of functions of L^(R+) with respect to two independent standard Brownian 
motions. This theorem permits us to obtain, in Section 3, results of convergence in 
law, in the space C([0,T]), towards fBm, the process defined in (|1.8p and, finally, 
sub-fBm with parameter H S (0,1) using the decomposition of this process as a 
sum of two independent processes. 

Positive constants, denoted by C, with possible subscripts indicating appropriate 
parameters, may vary from line to line. 

2. General convergence result 

In this section we prove the main result of weak convergence in the sense of 
the finite dimensional distributions. We will use this result later in order to prove 
weak convergence results towards fractional Brownian motion and sub-fractional 
Brownian motion. 

Theorem 2.1. Letf{t,-) andg{t,-) be functions of L'^{R+) for all t e [0,r] ,r>0, 
let {Ns,s > 0} be a standard Poisson process and 6 G (0,7r) U (tt, 27r). Define the 
processes and , which are given by Y-f — {J^^ f{t, s)dWs,t G [0,r]} and 
fs = {J^ g{t,s)dWs,t G [0,r]} and where W = {W,,s > 0} and W = {W,,s > 
0} are independent standard Brownian motions. We also define the following pro- 
cesses 

(2.1) F/ = |^^ fit, s) cos (ON^') ds, te[0,T] 
and 

(2.2) f/ = 1^ g{t,s) sin (^0N^^ ds, t <E [0,T] 

Then, the finite dimensional distributions of the processes {Y/} and {Y^} con- 
verge in law to the finite dimensional distributions of the processes Y^ and Y^ . 

Proof. Taking into account that the proof is valid for any fixed t G [0, T], by abuse 
of notation we will write /(s) instead of f(t,s). Slightly modifying the proof of 
Theorem 1 in |D JOO] . in order to prove the weak convergence, in the sense of the 
finite dimensional distributions, it suffices to show that 

(2.3) E[{Y/f] <c(^J^ f^{s)dsy E 



9^2 



< C 



g\s)ds 
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Observe that defining 

f f ~ f 2 f°° , , i0N2s 

Zi = Y/+iY/ = - f{s)e l^ds 

£ JO 

we have E[Z/Z/] = E[{Y/)^ + {Y/f]. Therefore if we prove < C||/||i, 

where ||-||2 is the L^(IR+) norm, the stated convergence follows. 



E[Zfzf] = E 



2 



ieN2 







/(s)e ^ds- / /(r)e ^ dr 



4^ 



OO /'OO 



Jo 

OO /"OO 



i0( N2s -N2t I 

f{s)f{r)e V T^/dsdr 



JO 

OO PQC 



+ 



Jo 



l{.<.}/(«)/WlE 
l{s<r}f{s)f{r)E 



dr ds 



01 N2r- -N2s 
■71 71 



dsdr. 



Since E[e*^^] = e-^^i-*'") and Ele"*^-^] = e-^^^"* being X a Poisson random 
variable of parameter A, we obtain 



A /'OO /"OO 
^ /"OO poo 

+ ^ / l{.<r}/(s)/(r-)e-2'?(^-^"' 'drds 
2^ Jo Jo 

<-/ / l{r<.}|/(s)/(r)|e-2^(i— ^^)drd5 
£ Jo Jo 

A POO POO 

+ -/ / l{.<r}|/(s)/(r)|e-''^^('-^°^'^drd5. 

^ Jo JO 

Using the inequahty \f{s)f{r)\ < \ {P{s) + P{r)) and noting that, by means of a 
change of variables, the last two integrals are the same we have that 

A POO POO 

nzfzf] < 4 / / {fis)+fir)) e-^^(i— ^)drds 

^ ^0 ^0 

A / poo poo poc pr N 

= 1// / e-^^(i— «%rds + / /2(r)/e-2^(i-™«^)dsdr 



< ^ 



1 — cos I 



1 



Then, the convergence of the finite dimensional distributions has been proved 
and it remains to prove the independence of the limit processes. We begin by 
proving that the family {Y/Y^}^^o is uniformly integrable. Indeed, we will prove 

that supg-^o 

< 00. Using Holder's inequality we have 

supE \{y/yM < sup (E[(y/)^])^ (E[(y/)^])' . 

£>0 L J £>o V / 
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In order to prove that the last expression is finite, we will show that 



(2.4) E[(r/f]<C7(^y^ f{s)dsj, E[(f/)4]<c(^y^ g\s) ds 
Being Z/ hke before, we can prove ([211) showing that E[(Z/Z/)2] < C\\f\ 



16, 



E[(Z/z/)2] = -E 



[0,00) 



i9 N2^ +N2^ -^2£, -A'2£4 

/(si)---/(s4)e V ^ ^ ^/dsi---ds4 



64 



/ l{si<...<.,}/(si) • • • /(S4)E [Si + • • • + £^6] dsi • • • ds4 

J[0,oo)4 



where 



( ^2 = 1 +N2S2 -N2S, 



-7^ -7^ 



E4 = E3, — E2, — El. To obtain the last expression note that we can 
arrange si, S2, S3, S4 in 24 different ways and due to the symmetry between si and 
S2 and between S3 and S4 we have 6 possible different situations, Ei, . . . , Eq, each 
one repeated 4 times. By means of the properties of Poisson process we have 



i[Ei]l\\E[E2]l\\nE3]\\< 



^-2 °^ /3 (i_cosg) -2 (l-cosfl) 



and we can conclude 

Eiizizir] < ^ / i{,,<...<,,}i/(si).../(s4)i 



< 



384 
2^ 



[0,00)2 



< 3 



1 — cos ( 



I{.,<..}l/(si)/(s2)|e-2'^(i-™««)dsids2 



Then the family {F/F/} £>o is uniformly integrable and consequently 



E 



lim r/(t)f/(s) ^ lim E[y/(t)y/(s)]. 



bmce y-'' and Y3 are centered Gaussian processes, in order to prove their indepen- 
dence it suffices to show that the last limit converges to zero as e tends to zero. To 
deal with this limit, we observe that 



00 />OQ 



2 , , /(s).g(r)E cos(07V2.)sin(^?7V2.) 
£ Jo Jo '- - ^ 



dsdr 



00 poo 



Jo 

00 POO 



/(s)5(r)l|,<,}E cos(0iV2.)sin(0A^2.) 



dsdr 



"'0 



/(s)5(r)l{,<4E cos{0N^)sm{eN2r) 



dsdr 



I! + 11 
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Applying the formula 2 sin a cos b ~ sin(a + b) + sin(a — 6) = sin(a + b) — sm{b — a) 
we have 

/•OO /'OO ^ 



^0 



/(s)ffWl{s<r}E sin(0(7V2. + )) + sin(0(iV2. - N..)) 



~ ^1,1 + -'1,2) 

/•OO /'OO 



"'0 
-'2,1 -'2,2 



/(s)g(r)l{,,<,}E sin(0(7V^ + A^- )) - sin(0(A,. - iV,.)) 



dsdr 



dsdr 



We proceed to show that l\ ^ and /| ^ converges to zero as e tends to zero and that 
1^2 and If 2 have the same (finite) limit, thus obtaining the stated result. We note 
that 

/"OC />oo 



I\. = Im 



/(s).9Wl{.<r}l 



iO{N2r -N2s) 2ieN2s 



dsdrj 



Im(A" 



To find the limit of ^ we sec that \\A^\\ converges to zero as e tends to zero. 

^•00 /"OO 

\-tl \ I M-f — %(r — s)(l — cos9) — it-sfl — cos 26) 1 j 

/(s)5('")|l{s<r}e ~ ^ Msdr 



i^l < ^ 



^0 

OO pOQ 



1 /"OO /"OO 

<J- / (/2(s)+g2(^))l^^^^je-fT(l-cose)3(cos29-cos0)j^^^ 

^ Jo Jo 



3 



r (cos 2^? — COS 6) 



r(l-cose) 



drds 



1 



.92(r)e-S(i- 



i-(cos2e-cose)j^^^ 



When cos 6* cos2(? it is easy to check the convergence to zero. Otherwise, we 
integrate obtaining 

1 



2(1 - cos 6*) Jq 



f-(l-cos2e) 



ds, 



AS — I 

2 2(cos26l-cos6i) 
1 



( 



(1-COS61) / ■25-(cos2e —COS 6 



/■OO 

y (e-i^(i-™^2e) _ g--(i-cos9)^ 



2(cos26' - cos 6*) 

which concludes, as the convergence to zero is easily seen by dominated convergence. 
In the same manner we can see that /f i converges to zero. 
With respect to the term If 2 we observe that 



00 />00 



^0 



/(s)5(Ol{s<.}e-*('^-^^(^-^"'^dsd7 



Since ^(1 — e**)e ^'■'^ "■'^^ is an approximation of the identity, we have that 
If 2 converges, as e tends to zero, to Im (^j^-iw fis)g{s)ds^ < 00. Clearly the 
same result is obtained for /| 2 • This finishes the proof. □ 

Remark 2.2. We can use this result to approximate two independent processes of 
many kinds, such as processes with a Gousart kernel (see for instance |DJOO| ). 
the Holmgren- Riemann-Liouville fractional integral ( |DJOO| ). fractional Brownian 
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motion and sub-fractional Brownian motion. In the next section we study weak 
convergence, in the space C{[0,T]), towards these two last processes. 

3. Weak approximation of fractional and sub-fractional Brownian 

MOTION 

In this section we apply the main theorem to prove weak convergence results 
towards fractional Brownian motion and sub-fBm. Let us begin with fBm and 
later on we will study the convergence for sub-fBm. 



3.1. Weak approximation of fractional Brownian motion. 

We are going to prove a result of weak convergence in C([0,T]) towards fBm, 
applying Theorem l2.1l In order to do so, we use the following representation of the 
fBm as the integral of a deterministic kernel with respect to standard Brownian 
motion (see for instance |DU99| ) 

(3.1) = f K"it,s)dWs, 



du, 



where H G (0, 2), {t, s) is defined on the set {0 < s < i} and is given by 
(3.2) K^it, s) = d"{t ^s)^+d^ (i^) - (l - (i) 



where the normalizing constant is 



T{S±1)T{2-H) 



2 

Since in this section the domain of fBm is restricted to the interval t G [0,T], we 
can rewrite the integral representation as 

Jo Jo 
where K"{t,s) = (t, s)l[Q,t^{s). 

Applying this representation, since {t, •) G L'^(M+), the following result is a 
corollary of Theorem 12.11 

Corollary 3.1. Let K"{t,s) = K" {t, s)l[o,t]{s), where K"{t,s) is defined by 
iS.S^) . let {Ns,s > 0} 6e a standard Poisson process and let 9 G (0,7r) U (7r,27r). 
Then the processes 

= 1^^ K" {t,s) cos (oN^'j ds, tG[0,r]| 



(3.3) 
and 

(3.4) I ^ K"{t, s) sin (ON^^ ds, t G [0, T] j 

converge in law, in the sense of the finite dimensional distributions, towards two 
independent fractional Brownian motions. 

We now proceed to prove the continuity and the tightness of the families of pro- 
cesses defined by p.3|) and ()3.4p . and consequently, proving the weak convergence 
in the space C{[0,T]). 

Theorem 3.2. Under the hypothesis of Corollary \ 3. 11 if moreover one of the fol- 
lowing conditions is satisfied: 
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(1) H€{l2), 

(2) H e (0, i] and satisfies cos((2i + 1)6*) / 1 for all i e N such that i < 
2 [ti] ' 

then the processes and converge in law in C([0, T]) towards two independent 
fractional Brownian motions. 

Proof. We first observe that the processes and B^ are continuous. In fact, B^ 
and B^ are continuous for all H G (0,2) and absolutely continuous if e (lj2), 
since it can be proved that (see Lemma 2.1 in [BF05j ) 

\B^(t)-B^{s)\<CH{t-si'^y^ 

and 

\Bf{t)-Bf{s)\<CH{t-s)i'^y\ 

It only remains to prove the tightness of the families of processes defined by (13. 3p 
and (|3.4p . Since B^ {0} = 0, using Billingsley's criterion (see for instance [Bil68| ) 
it is enough to check that for some m > and a > 1 

E[|5f (i) - B^{s)r] < C{F{t) - F{s)r, 

where -F is a nondecreasing continuous function. 
On the other hand, it is known that 

r {K^{t, r) - K"{s, r)ydr = E [(Bf - B^f] = {t - s)^ 
Jo 

and then it is sufficient to show that 

(3.5) E[{yfr]<Cra(^£ fir)dr^ , E [{yfr] < f\r) dr^ 

holds for some m satisfying the condition Hm/2 > 1, where /(r) :— K^{t,r) — 
K"is, r), yf = f fir) cos(0iV|. )dr and g = f f{r) sm{eN^)dr. 

Then, in the case ([T]), it is sufficient to prove (|3.5p for m — 4, which can be seen 
proving that E[(z/z/)^] < C||/||2, where ||-||2 is the i^[0, T] norm and z/ = y^+iy{. 
If we extend / to 1R+ for zeros, i.e., if we consider F{r) :— f{r)t[Q^T]{i^)j we have 
proved in Theorem 12.11 that 



E[(Zfzf)^]<3(^^%^(.)d. 



2 



Then, 

E[(z/z/)2]=E[(Zfzf)2] 



T \2 



To prove the result under the hypothesis ^ we can show that (|3.5p is satisfied 
for some even m such that > 1. If we proceed in the same way as in case ([T]) 
we obtain an expression that depends on 1 — cos((2i + 1)6*) for alH = 0, 1, . . . , [^] 
and the constant Cm depends on max^^g ^ i_cos((2i+i)e) • '-' 
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3.2. Convergence towards sub-fractional Brownian motion. 

In order to obtain the convergence to sub-fractional Brownian motion, we will 
apply a decomposition result due to Ruiz de Chavez and Tudor in [RdCTOQ] . In 
this paper, they use the process introduced by Lei and Nualart in [LNOQj and 
defined in (|1.8p by the equation 

/•oo 

= / (f -e-''*)r-^diy^, 
Jo 

where is a standard Brownian motion. It can be proved (see [LNOQj or [RdCTOO] ) 
that its covariance function is 

r3 6) Co.(X" X") / ^ ^'"^ + -s"-it + s)"] if H e (0, 1), 

(3.6) Cov(X, , X, ) - j ^ ^ ^ 

and that X^ has a version with absolutely continuous trajectories on [0, oo). 
The decomposition result can be stated and proved as follows: 

Theorem 3.3 (Decomposition of sub-fBm). Let be a fBm, a sub-fBm and 
W — {Wt,t > 0} a standard Brownian motion. Let X^ be the process given by 
m.8]) . If for H G (0, 1) we suppose that B^ and W are independents, then the 
processes {Y^ = CiX^ + B^ ,t > 0} and {S^ ^t > 0} have the same law, where 

^1 ~ \J 2T(i- g y • V for H G (1,2) we suppose that and W are independents, 
then the processes {Y^ — C2X^ + , t >Q} and {B^ , t > 0} have the same law, 
where C2 ^ "rS-H) ■ 

Proof. It is clear that the process Y^ is centered and Gaussian in both cases. For 
H G (0, 1), from (ll.6|l . ()3.6p and using the independence of X^ and B^ we have 

Cov(r,^, y/) = cfcovix^, ] + Cov[b^, sf ] 

= s" + t"-'-[is + t)" + \s-tn, 

which completes the proof in this case, and for H G (1,2), from p.7p . (|3.6p and 
using the independence of X^ and we have 

Cov(y,^, y/) = c|Cov[xf , ] + Cov[5f , s^] 

= l{s"+t"-\s-tn, 

which completes the proof. □ 

In order to apply the main theorem to prove weak convergence to sub-fBm, we 
have to prove weak convergence to fBm and the process X^ introduced by Lei and 
Nualart. Then, applying the decomposition theorem and the independence of the 
limit laws, we can state the weak convergence to sub-fBm for H G (0, 1). 

So, it just remains to prove for the process X^ defined by (jl.Sp the same results 
we have obtained for fBm. 

Corollary 3.4. Let X" be the process defined by iL8\) . let {Ns, s > 0} be a standard 
Poisson process and let 6 G (0, tt) U (tt, 27r). Then the processes 

(3.7) Xf = |^^ {1 ~e-'*)s-^ cos (eXgs^ ds, tG[0,r]| 
and 

(3.8) lf = |^^ {1 ~e-'*)s-^ sin (^0X2^'^ ds, iG[0,T]| 
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converge in law, in the sense of the finite dimensional distributions, towards two 
independent processes with the same law that . 

Theorem 3.5. Under the hypothesis of C'orollary \3.4\ the processes defined by Jg. ?| ) 
and i3.8]) converge in law, mC([0,T]), towards two independent processes with the 
same law that the process defined by il.8\) . 

Proof. We first need to show that the processes X^ and X^ are continuous. In 
fact, they are absolutely continuous. Let us consider for all r > the process 

2 



Yr. 



s 2 e COS 



ds. 



This integral exists because, using inequality (|2.3p . we have 



< C 



,l-H„-2sr, 



ds] ^C; 



r(2 - H). 



On the other hand. 



f \Yr\dr < f (E[i;2])^dr < C / r^dr 
Jo J Jo Jo 



< oo 



since H e (0,2). 

Let us now observe that X^ = J^^ Yrdr. Indeed, applying Fubini's theorem, 



Yrdr = 



- / (1 - e^"*) cos [ONg^j ds 

H 



The same proof shows that the process X^ is continuous. 

Next, we prove the convergence only for p.7p . For (j3.8p the result is proved 
similarly. 

It suffices to prove the tightness of the family {X^}^. Since X^ [0) = 0, using 
Billingsley's criterion we only need to prove that 

¥.[\X^{t)-X^{st]<\F{t)~F{s)\^ 

where _F is a continuous, non-decreasing function. We observe that 

E[|Xf (t)-Xf(s)|'*] =E -/ {^"{t,r) ^^"{s,r))cos{eN^)dr 

Jo ^ . 

where <^>"{t,r) = (1 - e~''*)r~^ e L2(M+). 

Since <f>^ £ L^(K+), applying the bound (12. 4p . which is proved in Theorem 12. 11 
we obtain 



[\X!'it)-X^is)\^] < C (^^" {<f"it,r) - $^(s,r))'dr^ 

= ^(/" ((1 - e-'''fr-^'+"^ + (1 - e-'-'^)2r-(i+^) 
-2(l-e-'^*)(l-e-'^^)r-(i+^))drJ . 



Using (|3.6p and assuming s < i we obtain for H G (0, 1) 
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E [|Xf (t) - (s)|4] < C {2{t + sf {2t)" - {2s)"Y 

<C{{2t)^ -{2s)")\ 

since s + t <2t. In the same way, if iJ e (1, 2), 

E [|Xf (i) - (s)|4] < C {{2tf + {2s)" - 2{t + sff 

<C{{2tr -{2s)")\ 

since s + t > 2s. In both cases we have proved the resuh with F{x) = {2x)" . □ 

FinaUy, we obtain the result of weak convergence to sub-fractional Brownian 
motion, as a direct conclusion of the previous results. 

Theorem 3.6. Let H G (0,1), let {X^{t),t e [0,r]} be the family of processes 
defined by (3 ^, let { i3f{t),t e [0,T]} be the family of processes defined by \3.4^ 
and Ci = \J 2r{i-il )' assume 9 S (0,7r) U (tt, 27r) and, for H e (0, 

that 9 is such that cos((2i + 1)9) ^ 1 for all i N such that i < ^ [-^j . Then, 
{Yf{t) = CiX^(t) + B^{t),t e [0,T]} weakly converges m C([0,T]) to a sub- 
fractional Brownian motion. 

Proof. Applying Theorems 13.21 and 13.51 we know that, respectively, the processes 
B^ and converge in law in C([0, T]) towards a fBm and the process defined by 
p.Sp . Moreover, applying Theorem 12.11 we know that the limit laws are indepen- 
dent. Hence, we are under the hypothesis of Theorem 13. 3) which proves the stated 
result. □ 

Remark 3.7. Obviously we can also obtain the same result using the families of 
processes defined by p.Sp and (|3.3p . 
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